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Abstract 

The ratio n/s, shear viscosity (77) to entropy density (s), reaches its local minimum at the 
(second order) phase transition temperature in a wide class of systems. It was suspected that 
this behavior might be universal. However, a counterexample is found in a system of two weakly 
self-interacting real scalar fields with one of them condensing at low temperatures while the 
other remains in the symmetric phase. There is no interaction between the two fields. The 
resulting n/s is monotonically decreasing in temperature despite the phase transition. 



I. INTRODUCTION 



What are the most perfect fluids in Nature with the smallest shear viscosity (rj) per 

entropy density (s)? Kovtun, Son, and Starinets (KSS) [1] suspected that they are a class 

of strongly interacting conformal field theories (CFTs) whose rj/s — 1/(4-71"). They even 

conjectured that l/(47r) is the minimum bound for rj/s for all physical systems. Ever 

since the KSS bound was proposed, much progress has been made in testing this bound 

fl fl 

and trying to identify the most perfect fluid (see [2j, [3j for recent reviews). It is found 
that rj/s can be as small as possible (but still positive) in a carefully engineered meson 
system , although the system is metastable. Also, in strongly interacting CFTs, the 
universal value rj/s = l/(47r) is obtained only in the limit of infinite N, with N the size of 
the gauge group, and infinite t'Hooft coupling limit |6j. 1/N corrections can be negative, 
however, [7|, |8| and can modify the rj/s bound slightly 

In the real world, the smallest rj/s known so far belongs to a system of hot and dense 
matter thought to be quark gluon plasma just above the phase transition temperature 
produced at RHIC Uj with rj/s = 0.1 ± O.l(theory) ± 0.08 (experiment ) 12j . A robust 
upper limit rj/s < 5x1/ (4ir) was extracted by another group 13J and a lattice computation 



ermi 



of gluon plasma yields rj/s = 0.134(33) |l4j. Progress has been made in cold unitary 
gases as well. An analysis of the damping of collective oscillations gives rj/s > 0.5 [15l 
ven smaller values of rj/s are indicated by recent data on the expansion of rotating clouds 



16|. 



17], ll8| but more careful analyses are needed [19|. 

Previous studies have given some clues about where to find the most perfect fluid in 
nature. The first one is to study strongly interacting systems because strong interaction 
generally implies small rj/s. The second clue can be found in a large class of systems where 
rj/s goes to a local minimum near the phase transition temperature (T c ) [20, |2l|, |25|. In 
particular, rj/s develops a cusp(jump) at T c for a second(first) order phase transition and 
a smooth local minimum for a cross over. This behavior is seen in QCD with zero baryon 
chemical potential 20J, [2JJ] and near the nuclear liquid-gas phase transition 22]. It is 
also seen in cold unitary fermi gases 2311. in H2O, N, and He and in all the matters with 



data available in the NIST database 



20 



22 



24J. Theoretically, these behaviors can be 



reproduced in controlled calculations of weakly interacting real scalar field theories 25]. 
Thus, it was speculated that this feature is universal. If this is indeed the case, then rj/s 
can be used to probe some parts of the systems which are hard to explore otherwise. For 



example, one can try to locate the critical point of QCD by measuring rj/s |22|, |26 |. 

In this paper, however, we present a counterexample of the rj/s behavior speculated 
above. In this model, rj/s does not go to a local minimum at the second order phase 
transition temperature. Our model is a mixture of two weakly self-interacting real scalar 
fields with one condensing at low temperatures while the other remains in the symmetric 
phase. There is no interaction between the two fields. The advantage of this model is 
that its 77/s can be computed reliably as in |25| because of the small couplings [27|. Other 
counterexamples have been asserted previously in literature. One of them is a a model 
calculation with a local minimum below T c 28]. In this model, large couplings are used 



to mimic the case of QCD. Thus, it is not clear this is due to the failure of the Boltzmann 
equation at large couplings 27|, or if the effect is generic. Also, holographic models 
have constant 77/s (= 1/4tt) in the limit of infinite N and the infinite 't Hooft coupling 
limit. If 1/AT corrections are added, rj/s becomes monotonically increasing below T c and 
a constant above T c [29] . Our model is a field theory model that one can compute directly 
and reliably. Our final result shows that rj/s does not have to develop a local minimum 
at T c . 



II. THE MODEL 



We will study real scalar theories in cases I- 1 1 1 : 



Ci — C\, 
Cu = £2? 

Cm = C x + C 2 , (1) 



^ = \{d^i? - (2) 



where 



The r]/s of cases I and II are well studied in [25] and we follow the treatment there. A 
and $ are renormalized quantities and the counterterm Lagrangian is not shown. The 
renormalization condition is that the counterterms do not change the particle mass and 
the four-point coupling at threshold. We will set < A <C 1 such that the systems are 
bounded from below and we can compute to leading order, 77 and s, in the Aj expansions. 
In case I, [i\ > and <f>\ stays in the symmetric phase. The resulting rj/s is monotonically 




FIG. 1: m q p vs. T for cases I (solid curve, without out symmetry breaking) and II (dashed 
curve, with symmetry breaking). Parameters can be in arbitrary units. 



decreasing in temperature (T) 25j. In case II, [x\ < 0. The 02 —02 symmetry 
is spontaneously broken below the phase transition temperature T c . The resulting rj/s 
is monotonically decreasing when T < T c and becomes monotonically increasing when 



T > T c . Also, rj/s forms a cusp at T c under the mean field approximation 251 ] . 

Because there is no interaction between <\>\ and 02 in case III, the entropy density is 
just the sum of the 0i and 02 entropy 

sin = sj + s H . (3) 

Analogously, in a linear response theory, the Kubo formula relates r\ to an ensemble 
average of a correlator 



ri = ~\f df f dt [ da; 3 ([T^(0),T^(x,t)]) , (I) 

^ J —OQ J GO J 

where T 2 - 7 is the spacial part of the off-diagonal energy momentum tensor. Tj J jj = Tj* +Tjj 
and ([3f (0),2g(x,t)]> = [<I?(0)> , (2§(x,0>] = 0, such that 

??777 = ??/ + r///. (5) 



The high T behavior of r)/s can be analyzed using the 1/T expansion as in Ref. 221 ] . 
By neglect the slow running of the coupling constant, the dimensionful quantities fif and 



T can only contribute to the dimensionless ratio rj/s through the pf/T 2 combination (note 
that it is p 2 , not pi that appears in the Lagrangian). As T — > oo, rj/s has the following 
1/T expansion 



,2 



? § l + c 2 ^ + 0(T- 3 )), (6) 



s/ A 2 V T 2 
^-§( 1+ 4 + (T-3,), (7) 

where ci > and c 2 > 0. There is no 1/T term because as mentioned above, the result 
does not depend on the sign of p^. The 1/T expansion of n/s in case III has the similar 
structure 

A 2 V 2 2T 2 ' 
Furthermore, in the limit of p\ = pf, we have rjj = n n , sj = Sji, and rjui/sjjj = rjj/si by 

Eqs.([n]) and §5§. This implies c[ = c\ and c' 2 = c 2 > 0. Therefore, if /z 2 + /if, > 0, r] in /s in 

is monotonically decreasing in T as T — >■ 00. The question is, whether this behavior 

persists from large T down to T c . Before answering this question numerically, we will try 

to understand the behaviors of r] and s separately. 

In Fig. 1, we show the typical m qPj i, the effective quasiparticle mass of 
function of T. When T = 0, m 2 I = n\ and rn 2 n = 2 |/x| | . Then m 2 pI increases for 
increasing T due to the positive thermal mass effect, while m 2 qp ji decreases to zero at T c , 
and then increases again at higher T. As T — > 00, m 2 qp I — rn 2 qp II = fif — /x| > 0. We have 
chosen the parameters such that m 2 p I > m 2 p n at all T, which gives sj < sn in Fig. 2. 
The cusp in s/j is barely visible. 

In Fig. 2, 77/T 3 is also shown. Its behavior is very similar to that of m qp . To further 
explore this relation, we use the kinetic theory approximation 

r] ~ pvl, (9) 

where p, v and / are the quasiparticle density, velocity and mean free path, respectively. 
Then using I ~ 1/nva, where n is the number density and a is the cross section between 
quasiparticles, we have 

— = (10) 
no a 

where e is the averaged quasiparticle energy. In a weakly coupled system, e can be 
approximated as a gas of free particles with mass m qp 

I d 3 p e qp f (e qp ) 



I d 3 P f (e qp ) 



[l + 0(A)] f (11) 
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FIG. 2: s/T 3 and 77/T 3 for cases I (solid curve, without out symmetry breaking) , II (dashed 
curve, with symmetry breaking), and III (dotted curve, the mixture of I and II). Parameters 
can be in arbitrary units. 

where e qp = ^Jp 2 + m 2 p and the Bose-Einstein distribution / (e qp ) = 1/ [e eqp ^ T — l). In 
two body collisions, a can be approximated as 



\ 2 

a ~ —5-. (12) 



Thus, 



6 3 
A eff 



(13) 



The effective coupling A e // is T dependent. The explicit expression for the scattering 



amplitude is [221 ] 



%T ~ 6A + (6A(0)) 2 



+ 5- + 



s-m\ v t-m\ p u- m? qp 



(14) 



When T ~ 0, s ~ 4m 2 , p and t ~ u ~ 0. However, t ~ u ~ causes no momentum 
redistribution and hence the t- and w-channels have no contribution to 77. Thus, we can 
approximate A e // as 

6A 2 (0) 2 

Ae//-A + 7 -^. (15) 

Under this approximation, A e // decreases smoothly from 2A at T = 0, to A at T c and 
stays constant above T c . Since A e // only varies by a factor 2, we can further approximate 
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FIG. 3: e 3 /(A 2 T 3 ) vs. T for cases I (solid curve, without out symmetry breaking) , II (dashed 
curve, with symmetry breaking), and III (dotted curve, the mixture of I and II). Parameters 
can be in arbitrary units. 

it by a constant A such that rj ~ e 3 /A 2 . As shown in Fig. 3, the T dependence of e 3 /A 2 is 
indeed qualitatively similar to that of rj in Fig. 2. 

Finally, we present the rj/s results in Fig. 4. They are qualitatively similar to rj/T 3 . As 
speculated above, r} in /sjjj is indeed monotonically decreasing both below and above T c . 
This is a counterexample to the previous speculation that 77/s goes to a local minimum 
at T c of a second order phase transition. 

There are some approximations that we have made in this calculation but none of them 
should change our conclusion qualitatively. The first one is the "Hartree approximation" 
that is used to neglect all the sunset diagrams below T c . This approximation is good 
when T ^> \ l l 2 T c . At lower T, the Hartree approximation is not reliable. However, as 
T — > 0, s approaches zero exponentially (the excitations are massive) while 7] approaches 
zero via power laws. As a result, t]/s is decreasing in cases I-III at low T. This feature 
is not affected despite the Hartree approximation used. The second approximation used 
is the mean field approximation. Unaccounted quantum fluctuations can make the result 
reliable in the region |T — T c \ /T c < O(A). However, this region can be made arbitrarily 
small by reducing A. 

Finally, the end point of a first order phase transition (called a critical point (CP)) is 
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FIG. 4: r]/s vs. T for cases I (solid curve, without out symmetry breaking) , II (dashed curve, 
with symmetry breaking), and III (dotted curve, the mixture of I and II). Parameters can be in 
arbitrary units. 

also a second order phase transition. This is a special kind of second order phase transition 
which can be modeled by an effective field theory whose <pf and <pf couplings vanishes at 
T c , so that the leading coupling is </>f . Hence the CP case is different from the cases we 
consider here. In some systems rj diverges weakly near a CP. 

III. CONCLUSION 

The ratio r]/s, the shear viscosity (77) to entropy density (s), reaches its local mini- 
mum at the (second order) phase transition temperature in a wide class of systems. It was 
suspected that this behavior might be universal. However, we have presented a counterex- 
ample made of a system of two weakly self-interacting real scalar fields with one of them 
condensing at low temperatures while the other remains in the symmetric phase. There 
is no interaction between the two fields. The resulting rj/s is monotonically decreasing in 
temperature despite the phase transition. 
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